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An  Analysis  of,  and  Some 
Observations  on,  Dynamic 
Fracture  in  an  Impact  Test 
Specimen 

Numerical  simulations  of  crack-propagation  histories  in  four  cases  of  dynamic  tear 
test  experiments  on  4340  steel  are  performed.  The  influence  of  the  loss  of  contact  of 
the  specimen  at  various  times  with  either  the  supports  or  the  tup  or  both  is  critically 
examined.  In  each  case,  the  variation  of  the  dynamic  K -factor,  for  the  simulated 
crack-propagation  history,  is  directly  computed.  The  results  are  discussed  in  the 
light  of  current  controversies  surrounding  the  dynamic  fracture  toughness 
properties  governing  crack-propagation  under  impact  loading.  It  is  concluded  that 
these  controversies  may  not  be  fully  warranted. 


ABSTRACT 

Numerical  simulations  of  crack-propagation  his¬ 
tories  in  four  cases  of  dynamic  tear  test  experiments 
on  4340  steel  are  performed.  The  influence  of  the 
loss  of  contact  of  the  specimen  at  various  times  with 
either  the  supports  or  the  tup  or  both  is  critically 
examined.  In  each  case,  the  variation  of  the  dynamic 
K-factor,  for  the  simulated  crack-propagation  history , 
is  directly  computed.  The  results  are  discussed  in 
the  light  of  current  controversies  surrounding  the 
dynamic  fracture  toughness  properties  governing  crack- 
propagation  under  impact  loading.  It  is  concluded 
that  these  controversies  may  not  be  fully  warranted. 

INTRODUCTION 

Until  recently,  for  stituations  governed  by 
email-scale  yielding,  it  was  thought  chat  the  govern¬ 
ing  criterion  for  elastodynamic  crack  propagation  un¬ 
der  Mode  I  plane  strain  conditions  can  be  written  as: 
Kj(vft)  »  Kiq(v)»  where  KI{)(v)  is  the  velocity-depen¬ 
dent  fracture  toughness  of  the  material,  which  was 
thought  to  be  a  "reasonable"  geometry- Independent 
material  property.  This  hypothesis  appeared  to  have 
been  validated  in  several  studies  related  to  dynamic 
crack-propagation  Initiated  under  quasi-static  Loading. 
In  the  analysis  of  such  cases,  both  "generation"  and 
"propagation"  calculations  were  employed:  in  the 
former  calculation,  the  experiment  ally  measured  crack- 
propagation  history  was  simulated  to  find  the  stress- 
intensity  factor  or  the  velocity-dependent  fracture 
toughness:  the  latter  calculation  was  used  in  either 

of  the  two-wavs  (1)  based  on  a  given  versus  v 
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relation  to  find  the  crack  propagation  history,  or 
(ii)  to  find  the  best  K^D  versus  v  relation,  the  cal¬ 
culated  crack-propagation  history  corresponding  to 
which,  agreed  beat  with  the  experiment.  The  remarkable 
success  of  these  calculations  appeared  to  indicate 
that  the  prediction  of  dynamic  crack-propagation  and 
possible  arrest  under  general  loading  conditions  may 
be  well  within  the  grasp  of  current  art  of  computa¬ 
tional  mechanics. 

Recently  some  work  has  appeared,  however, 
that  seemed  to  cast  doubt  on  the  concept  of  dynamic 
fracture  toughness  that  is  independent  of  the  rate 
of  applied  loading.  In  Ref.[l],  Kanninen  et  al  re-  * 
ported  experimental  and  numerical  results  on  dyanamic 
tear  test  specimens  of  4340  steel,  a  high-strength,- 
rate-insensitive  material.  In  these  experiments, 
crack-propagation  was  initiated  from  notches,  of  ' 
varying  degree  of  "bluntness" ,  under  Impact  as  well 
as  quasi-scatic  loading.  The  dynamic  f racture- tough¬ 
ness  was  attempted  to  be  inferred  from  the  energy 
measured  to  be  absorbed  in  the  tear  test.  A  series 
of  "propagation  type"  linear  olasto-dynamic  analyse*, 
using  hypothetical  K  values,  were  performed.  A 
surprising  finding  ot  [l|  was  that  the  dynamic  frac¬ 
ture  toughness  governing  crack-propagat ion  initiated 
from  a  blunted  crack-tip  under  Impact  loading  m.iv  be 
significantly  higher  (roughly  170)  than  when  crack- 
growth  Is  initiated  cjiiasl-statlcallv  (rough  Lv  about 

65MNm“U3). 

The  primary  objective  of  the  present  paper 
is  an  attempt  to  unnivze  the  data  presented  in  [Ll 
and  to  examine  the  results  in  the  light  of  the  con¬ 
clusions  presented  in  ( 1 1 .  Tn  addition  to  the  blunt- 
notch  specimen  data  reported  In  [i|,  analysis  is  per¬ 
formed  also  of  the  data  for  a  fatigue  prc-c racked 
specimen  supplied  to  the  authors  bv  the  Batteile 
Columbus  Laboratories  in  January  1981.  In  the  pres¬ 
ent  paper,  in  contrast  to  those  In  (l],  "generation' 
type  calculations  are  employed,  i.e.,  the  experimen¬ 
tal  data  cor  'rack-velocity  versus  time  historv  is 
simulated  In  a  t  tnlt  *  element  program  to  determin- 


directly  the  stress-intensity  factor  variation  with 
time.  The  employed  finite  element  method  is  the 
"moving-singularity1*  procedure  reported  earlier  by 
the  authors  [2,  3].  In  the  present  analysis,  care¬ 
ful  attention  is  paid  to  the  boundary  conditions  on 
the  specimen,  especially  the  loss  of  contact  of  the 
specimen  at  various  times  with  either  the  supports, 
or  the  tup,  or  both.  Four  different  cases  of  experi¬ 
mental  specimen  are  analyzed.  In  each  case,  the  In¬ 
put  energy,  kinetic  energy,  elastic  energy,  and  frac¬ 
ture  energy  variations  are  computed.  It  is  noted 
that  each  ni*  these  energy  quantities  is  computed  di- 
rect  ly  in  the  present  procedure.  That  these  energy 
quantities  "blance  out"  is  nothing  but  an  a  posteriori 
check  on  the  present  calculations.  It  should  thus  be 
noted  that,  in  contrast  to  the  “propagation"  calcu¬ 
lations  in  [1],  fracture  energy  is  not  inferred  from 
an  energy  balance,  but  directly  computed  in  the  pres¬ 
ent  “generation"  calculation. 

Detailed  results  are  presented  for  each  of  the 
four  cases  analyzed.  These  results  are  analyzed  to 
arrive  at  some  “plausible"  conclusions  which  appear 
to  be  at  variance  with  the  conclusions  presented  in 

HI. 

ANALYSIS 

The  test  specimen  geometry  is  Indicated  in  Fig. 

1,  along  with  the  finite  element  mesh  employed  in  the 
modeled  portion  of  the  specimen.  Points  L  and  S  in 
Fig.  1  represent,  respectively,  the  loading  and  sup¬ 
port  points.  Sixty- two  8-noded  isoparametric  ele¬ 
ments  and  one  moving  singularity  element  are  used. 

The  specimen  geometry  indicated  in  Fig.  1  corresponds 
to  that  reported  in  [1],  and  a  plane-strain  condition 
is  invoked  in  the  present  two-dimensional  analysis. 

In  simulating  the  experiments  (1),  the  following  in¬ 
itial  conditions  are  used  in  the  present  analysis: 
at  time  t*0,  velocity  u^«6. 88m/sec.  The  tup  displace¬ 
ments  are  calculated  by  u^-^t.  The  present  analysis 

does  not  account  for  the  elasticity  of  either  the  tup 
or  the  supports,  largely  due  to  the  lack  of  knowledge 
to  quantify  such. 

In  all  but  one  of  the  present  four  series  of 
calculations,  account  is  taken  of  the  possibility 
of  lack  of  contact  of  the  specimen  with  either  the 
tup  or  the  supports  (i.e.,  the  tup  and  supports  can 
“push"  the  specimen  but  not  “pull")  at  various  in¬ 
stants  of  time,  as  and  when  the  analysis  may  naturally 
dictate.  In  one  case,  to  study  the  effect  of  the 
above  contact/no-contact  conditions,  the  specimen  was 
held  “fixed"  (i.e.,  the  tup  and  supports  are  always 
in  contact  with  the  specimen). 

Also,  it  is  to  be  understood  chat  the  present 
series  of  computations  are  the  so-called  “simulation** 
or  “generation"  studies  in  the  sense  defined  iu  (l|, 
as  opposed  to  the  “propagat ion"  studies  performed  in 
(1|.  To  the  uninitiated,  these  terms  imply:  "gener¬ 
ation  study"  means  solving  lor  the  dynamic  stress- 
In tens  try  tactur  for  a  crack  which  is  "made"  to  propa¬ 
gate  with  the  experimenat ily  measured  crack-length 
(or  crack-velocity)  versus  time  history,  while  "propa¬ 
gation  study"  implies  soving  for  the  crack-length 
(or  velocity)  versus  time  history  using  a  “hypothe¬ 
sized"  dynamic  fracture  toughness  (which  may  or  may 
not  depend  on  c rack- vo loi: 1 1 y  (l))  value. 

The  present  series  of  computations  are  sum¬ 
marized  in  Table  l. 


Table  I:  “Generation  Studies** 


Data 

Study  1 

No 

Notch-root 

Diameter 

Initiation 

Time 

Boundary 

Conditions 

D  TT  l 

0.064mm 

95  user. 

Fixed 

DTT  2 

0.064mm 

95  usee. 

Contact/ 

No-contact 

DTT  3 

0.064mm 

35  usee. 

Contact/ 

No-contact 

DTT  4 

0.000mm 

92.24  usee. 

Contact/ 

No-contact 

The  crack-length  versus  time  histories  for 
the  above  four  cases  are  shown  In  Fig.  2.  For  the 
cases  DTT  1  and  DTT  2,  the  Aa  vs  t  curves  used  are 
the  same  as  that  In  Fig.  4  of  Ref.  [1],  except  that 
the  boundary  conditions  are  different  as  in  Table  1. 
In  Fig.  4  of  Ref.  [1],  one  date  point  indicating  a 
crack-growth  of  -  2.5mm  at  t-45  usee,  is  included, 
but  this  point  was  omitted  in  [1]  In  plotting  the 
Aa  vs  t  curve.  In  the  present  study  designated  as 
DTT  3  above,  this  data  point  was  Included  in  the 
Aa  vs  t  curve,  and  further,  the  time  of  initiation 
of  propagation  was  chosen  such  that  the  K  value  at 

initiation  -  65MNm  *  ;  this  curve  is  shown  in  Fig. 

2.  DTT  4  indicates  the  data  obtained  from  M.F. 
Kanninen  [4)  for  a  fatigue  pre-ctacked  specimen; 
however,  the  Initiation  time  was  determined  to  be 
92.24  sec.  by  extrapolation  of  the  supplied  (41 
experimental  data. 

Prior  to  the  presentation  of  the  results, 
we  Indicate  briefly  the  analysis  procedure.  As  noted 
earlier,  the  present  analysis  of  dynamic  crack  prop- 
atation  is  based  on  the  procedure  developed  by  the 
authors,  and  detailed  elsewhere  (2,  3).  To  supple¬ 
ment  the  mathematical  procedures  in  [2,  31  for  the 
present  case,  we  consider  some  details  of  imposing 
“contnct/no-contact"  boundary  conditions  on  the 
specimen . 

We  designate  the  force  with  which  either  the 
tup  or  the  supports  "push"  the  specimen  as  being 
(+  ve).  Using  the  standard  notation,  the  reaction 
forces  at  the  points  where  displacements  are  pre¬ 
scribed  are  calculated  by: 

P  -  Kq  +  mq.  (1) 

-  a  '  *  ~ 

The  displacement  u  and  reaction  force  P  in  the  time 
step  (n+l)  ’are  predicted  by: 


Vl 


+  At 


n+l  n 


(2) 


(Pn‘Pn 

and  P  . ,  -  P  +  — L-  At  . 

n+L  n  At  n+1 

n 


O) 


It  is  noted  that  we  may  use  At  .  *At  *At.  Assume 
7  n+l  n 

that  P  ,  and  P  are  (+  ve)  and  that  P  f>P^.  Tbe 
n-L  n  n-1  n 

no-contact  condition  during  the  time  Increment  (n)  * 
(n+L)  is  predicted  to  occur  after  the  sub-increment 
of  time: 


If  0<At  <At,  during  the  (n+1)  step,  we  change  At  to 
c 

Atc  and  perform  the  analysis  with  condition  of  contact 
and  during  (n+2)  step,  we  change  At  to  Atp  (Atc+Atp 

•At),  and  perform  the  analysis  with  the  condition  of 
no-contact.  This  process  is  repeated. 

An  analogous  scheme  is  used  to  predict  the  tran¬ 
sition  from  a  "no-contact"  to  "contact"  condition; 
however,  this  time  by  monitoring  the  displacements 
of  the  respective  points  of  the  specimen  relative  to 
either  the  supports  or  the  tup. 


RESULTS 


First  we  consider  the  DTT  1  case.  In  this  case, 
the  displacement,  velocity,  and  acceleration  at  the 
point  L  (see  Fig.  1)  are  prescribed.  L^,  and 

as  shown  in  Fig.  3  are  the  times  when  the  reaction 
force  at  the  tup  becomes  zero  (note  "+  ve  reaction" 
implies  that  the  tup  is  pushing  the  specimen).  Nega¬ 
tive  "reaction  force"  is  observed  during  times 
Lj< t< ,  and  t>L^.  This  phenomenon  can  also  be  ob¬ 
served  in  the  experiments  of  Mall  et  al  [5].  Since 
the  tup  contact  time  measured  in  the  experiment  [1] 
was  about  180  usee.,  the  computation  was  stopped  at 
around  this  time.  Fig.  4  shows  the  variation  of  the 
computed  stress  intensity  factors  in  the  present 
"generation"  or  "simulation"  study.  The  times  marked 
by  fv^'  and  ^2'  In  Fig.  4  are  the  times  when  the 

crack  propagates  with  constant  velocities  v^  and  V2 

respectively  (see  [1]).  The  initiation  toughness 
(the  terminology  is  summarized  later  in  this 


paper)  obtained  in  this  computation,  as  seen  from 

Fig.  4  is  about  106MNm  *  .  The  computed  variation 
—  of  input,  strain,  kinetic,  and  fracture  energies 
with  time,  are  shown  in  Fig.  5.  It  is  noted  that 
In  the  present  procedure  [2,  3)  the  dynamic  K- factors 
are  solved  for,  directly.  From  this,  the  energy  re- 
lease-rate  Is  calculated.  Alternatively,  fracture 
energy  is  also  calculated  directly  fr or  a  crack-tip 
integral  of  work  done  in  separation  of  crack-faces. 
These  two  procedures,  discussed  in  (61,  were  noted 
to  give  almost  identical  results  for  the  present 


cases).  It  Is  noted  that  each  of  the  four  energies. 


Input,  strain,  kinetic,  and  fracture,  are  calculated 
Independently,  in  the  DTT  1  case  as  well  as  the  other 
three  cases.  In  the  present  work.  That  these  energies 
"balance"  is  an  Independent,  a  posteriori  check  on 
the  present  calculations.  It  is  seen  from  FLg.  5 
that  during  the  periods  L^<t<L2  and  t^L^  input  energy 


appears  to  actually  decrease,  due  to  "negative  reac¬ 
tion  forces"  during  these  times  as  discussed  earlier. 
We  will  comment  further  on  those  energy  variations 
later.  Fig.  6  shows  the  crack  opening  profiles  at 
various  times.  As  seen,  the  profiles  arc  nearly 
linear  except  very  near  the  crack-tip.  This  suggests 
the  possibility  of  developing  K  measurement  tech¬ 
niques  from  crack-opening  displacements,  in  dynamic 
tear  testing.  Fig.  7  shows  the  contours  of  equivalent 
stress  In  the  presently  used  singular  element. 


which  is  shown  hatched  In  FLg.  1.  As  seen,  the  stress 
level  reaches  that  of  yield,  j  ,  only  within  the  con¬ 
tour  shown  by  a  broken  line  In  Pig.  7.  The  maximum 
distance  of  this  cotour  Is  roughly  0.5mm  from  the 


crack-tip;  thus  indicating  the  insignificant  role  of 
plasticity  in  the  present  problem.  Contours  of 
principal-stress  difference  (which  may  be  used  in  a 
Tresca-type  yield  condition)  shown  in  Fig.  8,  can  be 
observed  to  be  more  or  less  similar  to  those  of 
equivalent  stress  (which  may  be  used  in  3  Mises-ty^e 
yield  condition)  shown  in  Fig.  7. 

We  now  consider  the  case  labeled  DTT  2  in 
Table  1.  It  is  noted  that  the  condition  of  contact/ 
no-contact  was  invoked  in  this  case.  As  seen  from 
Fig.  9,  the  specimen  is  not  in  contact  with  the  tup 
during  the  periods  L^<t<L25  L^<t<L^,  and  t>L^  as 

marked.  Also,  it  can  be  seen  from  Fig.  10  that  the 
specimen  is  not  in  contact  with  the  supports  during 
the  times  Sj<t<S2*  and  t>S^  as  marked  in  Fig.  10. 

Comparing  Figs.  9  and  10  it  is  seen  that  the  maximum 
reaction  force  Ps ,  at  the  support  point,  is  very 

close  to  the  maximum  tup  load,  P^.  The  variation  of 

the  computed  dynamic  K-factor  is  shown  in  Fig.  11. 

It  is  seen  that  the  initiation  toughness,  K_,  is 
-1  5  Id 

again  about  108MNm  *  .  However,  prior  to  initiation, 

K_  value  appears  to  reach  122MNm  *’^(>K  )  at  t*82 

1  id 

Usee.  Comparing  Figs.  9  and  10,  it  is  seen  that 
during  the  times  L^<t<S2*  S^<t<L^,  and  t>L^,  the 

specimen  is  not  in  contact  with  either  the  tup  or  the 
support,  i.e.,  the  specimen  is  a  free-flying  object  I 
Fig.  12  shows  the  variation  of  the  four  energy 
quantitites;  input,  kinetic,  elastic,  and  fracture. 
During  L^<t<L2»  L^<t<L^  and  t>L^,  since  the  specimen 

looses  contact  with  the  tup,  no  increase  in  input 
energy  occurs,  as  seen  from  Fig.  12.  It  is  noted 
from  Fig.  12  that  the  total  work  done  at  »180  usee, 
was  about  53  Joules.  This  is  less  than  half  of  the 
experimentally  measured  absorbed-energy  value  of  130 
Joules.  However,  in  the  present  analysis,  no  account 
is  taken  of  energy  dissipated,  if  any,  in  the  supports 
or  the  tup.  This  discrepancy,  while  of  no  consequence 
in  the  present  "generation-type"  study  (wherein  energy- 
balance  calculations  are  "by-products"  of  the  analy¬ 
sis),  can  have  deleterious  consequences  on  a  "prop- 
atation"-type  study  as  in  [1]  wherein  the  experimen¬ 
tally  measured  absorbed-energy  is  used  in  inferring 
the  fracture  energy,  from  which  a  "plausible"  frac¬ 
ture-toughness  value  is  calculated  [11.  This  is  an 
important  point  to  remember  in  comparing  the  present 
results  and  those  in  (1J.  In  Fig.  13,  the  crack¬ 
opening  displacements  in  cases  DTT  1  and  DTT  2  are 
compared.  It  i3  seen  that  the  COD  in  DTT  2  oscil¬ 
lates  around  the  values  in  DTT  l  case. 

Now,  we  consider  the  DTT  3  case.  The  tup 
contact-force  variation  is  shown  in  Fig,  14.  Com¬ 
paring  Figs.  9  and  14,  it  is  seen  chat  during 
L^<t<L?,  the  separation  between  the  specimen  and  the 

tup,  bigger  than  in  the  DTT  2  case.  This 

is  attributed  to  the  higher  compliance  of  the  DTT  3 
specimen  due  to  the  fact  that  growth  initiation  oc¬ 
curs  much  earlier.  Also  the  second  loss  of  contact 

of  the  specimen  and  the  tup  (L»«.t<L,)  and  the  second 

1  "» 

peak  of  P  (L  <t<L  )  are  smaller  than  those  in  the 
L  _  J 

case  of  the  DTT  2  specimen.  Fig.  15  shows  the  dis¬ 
placement  of  the  specimen  from  the  support  noint, 
and  the  support  reaction  force.  Comparing  Ft  4s.  10 

and  L5  it  is  seen  that  the  periods  S.-S,  and 

12  2  ■> 

are  longer  than  those  of  the  DTT  2  specimen;  ind  the 
peak  value  of  P^  is  smaller  than  that  of  DTT  2.  Again, 


3 


r 


these  tendencies  can  be  attributed  to  the  earlier 
crack  initiation  in  this  specimen.  Fig.  16  shows  the 
K-factor  variation.  Note  that,  as  shown  in  Table  I, 
the  crack  growth  initiation  time  was  chosen  to  be 

t*35  usee,  such  that  K  *65MNm  In  spite  of  this, 

la 

the  K- factor  variation  in  DTT  3  is  more  or  less  iden¬ 
tical  to  that  in  DTT  2  (Fig.  11)  until  about  t-95usee. 
During  the  period  of  i  95  usee,  to  146  usee.,  the 
K  value  in  DTT  3  decreases  while  that  of  DTT  2  in¬ 
creases.  The  variation  of  total  work,  strain  energy, 
and  fr?  ure  energy  are  shown  in  Fig.  17.  It  is  seen 
that  during  v^<t<v^»  a  very  small  amount  of  energy 

is  consumed  in  the  fracture  process. 

Finally,  we  consider  the  DTT  4  case.  The  tup 
contact-force  variation  is  shown  in  Fig.  18.  As  seen 
from  Fig.  2,  the  crack  propagates  faster  in  DTT  4 
case  than  in  the  DTT  2  case.  Due  to  this  reason,  the 
changed  compliances,  the  second  peak  of  (L^tcL^) 

is  smaller,  and  the  second  separation  between  the 
specimen  and  the  tup,  i.e.,  (u^-u^)  be tween ( L^< t<L^) , 

is  larger  in  the  DTT  4  case  than  in  DTT  2 .  The  sup¬ 
port  reactions,  and  separation  between  supports  and 
specimen,  are  shown  in  Fig.  19.  These  are  more  or 
less  similar  to  those  in  the  DTT  2  specimen.  From 
the  K-factor  variation  shown  in  Fig.  20,  it  is  seen 

that  the  K_.  value  is  about  lllMNrn  After 

la 

initiation,  K.  drops  significantly.  The  energy 
variation  plots  are  given  in  Fig.  21,  from  which  it 
is  seen  that  the  total  energy  to  the  specimen  is 
lower  than  in  the  other  three  cases. 

In  the  following  we  attempt  to  draw  some  rational 
conclusions  from  the  above  presented  numerical  data. 


Comments  on  Fracture  In  Impact  Specimens 

At  the  outset,  we  define  the  following  nomen¬ 
clature: 


K^:  Plane-strain  fracture  toughness  under  quasi- 

static  loading 

Kjjji  (Variable)  Dynamic  fracture  toughness  for 

a  propagating  crack 

K_d:  Initiation  fracture  toughness  under  dynamic 

loading 


kQ 

Qd 


"Apparent**  fracture  toughness  at  initiation 
with  a  blunt  notch  under  quasi-static 
loading 

"Apparent"  fracture  toughness  at  Initiation 
with  a  blunt  notch  under  dynamic  loading 


KID,min:  Louer  bound-  if  any>  for  KtD‘ 


A  typical  K  versus  c rack-veloc Ity  curve  widely 
reported  in  literature  is  given  in  Fig.  22.  For  must 

brittle  materials,  Kf_  appears  to  be  lower  than 
te  a.  ID, min  1  1 


K 


ID, min 


0.5  to  1,0  Klc. 


<S> 


The  governing  equation  tor  e las tody namic  crack  propa¬ 
gation,  currently  considered  to  be  valid  under  quasi- 
static  loading  conditions,  can  be  expressed  bv: 


KjU.v)  -  Kid(v) 


(6) 


where  K^D  is  a  function  of  crack-velocity,  in  general, 

the  above  equation  la  inconsistent,  from  a  functional 
view-point,  at  the  point  of  crack-growth  initiation. 
For  example,  if  the  "crack"  starts  propagating  (i.e., 
attains  a  finite  velocity  in  zero-time)  from  a  blunt 
notch  with  an  intensity  value  K^,  we  have: 


K 


I 


KQ(b)  +  K 


ID,min 


where  b  is  the  notch- root  diameter.  Moreover,  for 
initiation  under  quasi-static  loading  from  a  sharp 
notch  or  a  fatigue  crack  (b*0),  wherein  the  crack 
is  assumed  to  attain  a  finite  velocity  in  zero  time, 
we  have: 


Kt  -  K  *1  K  . 

I  IC  ID,  min 


(8) 


except  in  the  special  when  KID  min"Kxc* 

Now,  under  conditions  of  impact  loading, 

the  initiation  fracture  toughness  K  ,  as  reported 

Id 

in  most  of  the  literature,  is  generally  lower  than 
Klc«  The  values  of  K^/K^  as  reported  in  literature 

are  summarized  in  Table  2. 


Table  2:  Initiation  Toughness 
_ Under  Dynamic  Loading 


Investigators 

Material 

Temp  °C 

L_ 

Load  Rate 
kI[MNm"1,5/ 
sec.  ] 

Kid/Kic; 

Shabbits  [7] 

A533B 

10 

104 

-0.67 

Kalthof f 
et  al  [8] 

Araldite  B 

Ambient 

104 

*0.73 

I 

Mall  et  al 
(5] 

i  Poly- 
carbonate 

Ambient 

6xl03 

'0.65 

Ireland  [91 

A533 

10 

106 

>0.63 

Ireland  [9) 

4340  Steel 

i  . 

20 

Pre-cracked 

Charpy 

-1.00 

Influence  of  the  rate  of  loading  on  is  shown  in 

Fig.  23,  which  is  taken  from  Refs.  (7,  9].  As  seen 
from  Fig.  23,  for  most  of  the  materials,  de¬ 
creases  with  increasing  K^.  Fig.  24,  taken  from  Ref. 

[10]  shows  the  comparison  of  K  ,  K  (subscript  I 

d  0 

omitted  due  to  the  fact  that  plane-strain  conditions 
were  not  validly  mot),  and  data  for  A533B  steel. 

From  table  2  and  Figs.  23  and  24,  one  may  observe 
that: 


Ku  *  Ktd(V  M  °-5  ’  u0  Ktc-  (9) 


Once  again,  Kq.  (6)  is  funct ional 1 v  inconsistent  at 
growth-initiation  associated  with  an  impact  loading. 
Furchrc,  tor  growth- in  it  iaion  r'rom  a  blunt  notch 
under  impact  loading,  we  have: 


Wb,(VKId  JmJ  KQd  *  KID.nl«Ul0' 


However,  utter  a  certain  amount  of  crack  propagation. 
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h)  nay  be  applicable,  in  .t  rate- insensitive  ma- 
•.-rial,  even  under  impact  loading,  although  Eq.  (6) 

\  expressed  for  quasi-static  loading.  Contrary  to 
*w£St  for  a  rate-sensitive  material  under  impact 
leading,  the  governing  equation  for  dynamic  crack 
V.oagation  may  differ  from  Eq.  (6),  and  may  be  re- 
.ritton,  possibly,  as 

KjU.v)  •  Kid(v,c)  or  Kid(v,K)  (11) 

vnere  i  is  the  strain-rate  in  the  vicinity  of  the 
crack-tip  (see  Fig.  25).  If  £  (or  k)  is  small,  or 
the  material  is  almost  rate-insensitive,  the  situations 
in  Fig.  25  a  and  b  become  practically  the  same. 

It  should  be  noted  that  all  the  analyses  in  Ref. 
[1]  are  of  the  "propagation"  type  in  the  sense  de¬ 
fined  earlier.  In  [1],  the  K-  versus  v  curve  for 
4340  steel  is  first  assumed  asT 

Kid  -  65  +  0.044  v  (12) 

where  K.^  is  in  MNm  and  v  is  in  m/sec.  In  [lj, 

the  range  of  crack-velocity  for  which  Eq.  (12)  is 
valid  is  not  indicated.  However,  in  general,  there 
should  be  an  upper-limit  value  for  velocity  as  shown 
in  Fig.  22.  Additional  evidence  of  an  upper  limiting 
value  for  velocity  is  shown  in  Fig.  26,  taken  from 

Ref.  [10],  Froa  Eq.  (12),  K__  .  -  65MNm-1,5  while 

id, min 

KIC  aPPears  t0  about  50  from  the  Damage  Tolerant 

Design  Handbook  [  11 },  and  Fig.  26.  Considering  the 
near  rate-insensiei ve  behavior  of  4340  steel,  it  then 

appears  that  Kld=  K IC*  KID> nln=  5°~  “HNm" 1 ' 5  • 

In  order  to  sustain  a  running  cracfe,  the  sur¬ 
rounding  elastic  field  must  produce  plastic  strains 
continually  near  the  advancing  crack-dp  adequate 
for  the  separational  process  (region  R  of  Fig.  25). 

If  KId  is  considerably  lower  than  and  the 

elastic  field  near  the  crack-tip  cannot  continually 
supply  energy  to  the  process  zone,  the  crack  will  be 
arrested.  In  the  above  case  if,  on  the  other  hand, 
the  elastic  field  can  continually  supply  energy  to 
Che  process  zone,  the  crack  will  propagate  with  a 
somewhat  low  velocity  in  the  earlier  stages  of  crack- 
propagation.  This  appears  to  be  verified  in  the 
present  set  of  simulations.  Comparing  Che  varia¬ 
tions  in  DTT  2  and  DTT  3  specimens  (Figs.  11  and  16 
respectively)  It  is  seen  that  the  variations  are 

essentially  similar  evenchough  the  crack  Ln  the  DTT 
3  case  has  been  assumed  to  propagate  with  an  inltla- 
*°l  .5 

tion  value  of  K  *65MNm  .  This  can  bo  attributed 
to  the  nature  or  the  presently  considered  experiment 
in  which  the  elastic  field  near  the  crack  Is  always 
growing  Ln  intensity  during  the  period  t-10  to  80 
usee. 

However,  in  the  DTT  4  specimen,  the  present 
simulation  indicates  initiation  thoughness  of 

-  ilOMNm  (see  Fig.  20).  However,  .is  mentioned 
earlier,  the  in  It  i  u  ion  time  for  the  DTT  4  case  which 
was  unknown  in  the  experiment,  was  determined  by 
extrapoi jtlon  of  the  crack-growth  versus  time  data. 
Thus,  as  .in  be  rationalized  from  a  comparison  of 
the  HIT  2  and  HIT  )  cases,  one  m.iv  surmise  that  crack- 
growth  Initiation  n.  iv  have  occur.d  much  bet  ore 
y*ec.  (as  used  in  the  present  cnlciil.u  lens,  see  Fig. 
20).  even  It  with  i  verv  smalt  velocity.  Hier»  tore, 
one  mav  postulate  chit  even  lor  the  DTT  «  case,  K 


value  may  actually  have  been  Jess  than  or  equal  to 

65MNm  This  hypothesis  may  in  fact  be  considered 

to  be  also  supported  by  the  resuLts  of  (11.  Refferlng 
to  Fig.  8  of  (1),  it  is  seen  Cha  in  the  "propagation" 
calculation,  the  use  of  the  relation  Kjp*654-0 ,044|v 

actually  gives  a  much  better  agreement  with  the  experi¬ 
ment  in  terms  of  crack- Initiation  time  rather  than  the 

use  of  Kj^-l70MNm  Note  that  the  experimental 

data  for  the  first  measurement  of  crack  grwoth  is 
missing  in  Fig.  8  in  Ref.  (1(.  Contrary  to  this, 
after  a  certain  amount  of  crack-growth,  use  of 

K.jp*170MNm  appears  to  give  a  good  agreement  of 

the  "propagation"  calculation  [1J  with  the  experimen¬ 
tal  data. 

However,  as  noted  earlier,  the  value  of  KT 
-1  5 

»170MNm  *  was  "derived"  in  (1)  from  an  energy  bal¬ 

ance  consideration:  the  total  energy  absorbed  during 
the  impact  fracture  of  specimen  was  130  Joules;  the 
energy  imparted  to  an  "intact"  unconstrained  specimen 
from  an  elastic  collision  between  the  top  and  speci¬ 
men  was  79  Joules;  and  the  difference  (130-79)-51 
Joules  was  assumed  to  be  consumed  as  fracture  energy. 
From  this  assumed  fracture  energy  and  average  frac¬ 
ture  toughness  value  of  170MNm  was  derived  in 
(1 ] •  Thus,  the  propagation  analysis  {!]  represents 
an  a  priori  "energy-balance"  condition,  in  which  other 
sources  of  energy  dissipation  are  ignored.  On  the 
other  hand,  as  already  noted,  in  the  present  "genera¬ 
tion"  calculations,  an  energy  balance  relation  is 
an  a  posteriori  "by-product"  of  the  calculation  it¬ 
self  in  which  the  experimental  crack-growth  history 
Is  simulated.  In  connection  with  Fig.  12,  (DTT  2 
specimen,  which  is  also  the  basis  for  the  computation 
K^*170  in  Ref.  [1])  it  is  noted  that  the  computed 

total  input  to  the  specimen  is  53  Joules,  which  is 
less  than  half  of  the  experimentally  measured  value 
of  130  Joules.’  This  brings  to  question  the  neglection 
of  ocher  dissipated  energies,  and  hence  the  value 

of  Kjp-l70MNm  * ’ ^  hypothesized  in  [l]. 

The  present  analysis,  and  the  above  dis¬ 
cussion,  appear  Co  lead  us  to  believe  that  the  use  of 

versus  v  curve  with  .  *  65MNm  and  K_.< 

ID  ^  -  ID, min  Id- 

65MNm  ’  #  with  the  curve  being  such  that  the  upper 
value  of  v  is  limited  (as  in  Fig.  22),  may  be  war¬ 
ranted.  Suppose  this  "saturated"  curve  is  used  tn 
the  "propagation"  calculation;  then  the  crack  will 
start  to  propagate  ac  t  *  25  *  35  usee.  During 
most  of  the  period  from  initiation  to  about  t  »  100 
usee.,  the  crack  will  propagate  with  a  relatively  slow 
velocity,  as  discussed  earlier  (l.e.,  crack  growth 
will  be  small  until  about  t  -  100  usee.).  The  crack 
will  attain  a  highcr  velocity  because  of  the  nature 
of  the  testing  system  considered,  in  that  is  in¬ 
creased  continually  by  the  dropping  tup.  Then  the 
crack  speed  will  saturate.  In  addition  to  the  prop¬ 
erty  of  the  "saturated"  K^vs  v  curve,  there  may  be 

other  f.ictors  to  limit  the  crack  speed  In  the  present 
type  of  specimen.  One  of  these  is  the  presence  of  trie 
compressive  stress  field  In  the  ligament  of  the  speci¬ 
men.  is  can  be  seen  from  the  moment un  ha  Lance  condi¬ 
tion. 

Cootrarv  to  the  impact  loading^  under  qutsi- 
st.itic  loading  the  analysis  using  ♦  J.  ^4  v 

gave  eood  igrcemunt  with  the  expedient  in  (ll.  In 
this  case,  *ince  the  tup  was  fixed  to  correspond  to 
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K^-10SMNm“1’5,  the  vaLue  of  K  ^  during  crack-propa¬ 
gation  under  quasi-static  loading  condition  is  always 
lower  than  108MNra”L*5  [Kj (t , v)<K^] .  This  suggests 
that  K^*65  +  0.044  v  is  valid  for  lower  values  of 


velocity,  while  this  linear  relation  may  be  invalid 
for  Larger  velocities,  and  eventually  the  velocity 
will  be  limited,  leading  to  a  satuarated  K^D  versus 
v  curve. 
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